Speculations on primordial magnetic helicity 
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Abstract 

We speculate that above or just below the electroweak phase transition mag- 
netic fields are generated which have a net helicity (otherwise said, a Chern- 
Simons term) of order of magnitude Nb + Nl, where Nb,l is the baryon or 
lepton number today. (To be more precise requires much more knowledge of 
B,L-generating mechanisms than we currently have.) Electromagnetic helic- 



es \ ity generation is associated (indirectly) with the generation of electroweak 



Chern-Simons number through B+L anomalies. This helicity, which in the 



early universe is some 30 orders of magnitude greater than what would be 
expected from fluctuations alone in the absence of B+L violation, should be 
reasonably well-conserved through the evolution of the universe to around 
the times of matter dominance and decoupling, because the early universe is 
an excellent conductor. Possible consequences include early structure forma- 
tion; macroscopic manifestations of CP violation in the cosmic magnetic field 
(measurable at least in principle, if not in practice); and an inverse-cascade 
dynamo mechanism in which magnetic fields and helicity are unstable to 
rS \ transfer to larger and larger spatial scales. We give a quasi-linear treatment 

of the general-relativistic MHD inverse cascade instability, finding substantial 
growth for helicity of the assumed magnitude out to scales ~ hie' 1 , where e 
is roughly the B+L to photon ratio and Im is the magnetic correlation length. 
We also elaborate further on an earlier proposal of the author for generation 
of magnetic fields above the EW phase transition. 
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I. INTRODUCTION 



The problem of the generation of cosmic magnetic fields is still open. 
Some authors |]J argue that no primordial seed field is necessary, and that 
Biermann-battery effects (non- vanishing Vn x VT where n is density and T 
is temperature) act as a source for magnetic fields at the time of structure 
formation, with the source due to shocks. But many others search for primor- 
dial seed fields, associated with the electroweak phase transition p], || |], [|], 
the QCD phase transition |], 0, or various other mechanisms, e.g. inflation 
||. Efforts involving phase transitions are further subdivided according to 
their assumption of a first-order phase transition [[| or otherwise; a first-order 
transition produces bubbles and turbulence. 

So far there has been little discussion of the generation of primordial 
magnetic helicity Hm, defined by what is otherwise known as a Chern-Simons 
term: 

E M = [ d 3 xA-B. (1) 



This P- and CP-odd function is important for several reasons: 

• It is nearly conserved in the early universe (exactly so if the conduc- 
tivity is infinite). 

• It is not possible to have magnetic fields which are completely homoge- 
neous and carry helicity, so if there is primordial helicity there is also 
some sort of spatial structure in the universe at early times. 

• If it were possible to measure a net helicity for the present universe it 
would be another macroscopic manifestation of CP violation.^] 

• It is well-known to students of magnetohydrodynamics (MHD) that the 
presence of magnetic helicity (or other parity-odd expectation values) 
can lead to unstable dynamo action [10, [LI, 0, [TT~ 



In the present paper we consider a scenario involving generation of mag- 
netic fields and helicity at the electroweak phase transition. We can dis- 
cuss the part involving just the generation of magnetic fields fairly precisely 

lr The most important manifestation is the very existence of the stars and galaxies at 
present abundances. Much more speculatively, another could be the predominance of life 
forms of a single chirality. 
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[0, [14], it involves a condensate of EW magnetic fields in the magnetic 
N = Matsubara sector of the SU (2) part of the EW gauge theory. This sec- 
tor of the theory is strongly-coupled and an SU (2) magnetic condensate can 
be shown, without any approximation, to form This condensate is char- 
acterized by a finite density of vortices (closed string-like objects which are 
very long, with a thickness inversely proportional to the dynamically-induced 
gauge-boson mass) which are randomly linked; the Chern-Simons number is 
a direct measure of this linkage |15|, [14|, [16], [17]]. (Thermal sphalerons may also 
be present, but their Boltzmann factor is likely to be small |18| .) Like any 
other condensate, the vortex fields are sustained by a macroscopic number of 
phase-coherent W-bosons. When the temperature falls substantially below 
the critical temperature T c , the condensate loses phase coherence and is ex- 
pressed as particles, plus Maxwell magnetic fields generated by the charged 
Ws. These interact electromagnetically, retaining their original structure for 
a short while (because the W-boson mass does not change very much in the 
immediate neighborhood of the phase transition). For these Maxwell fields 
the Chern-Simons number or helicity expresses a linkage between magnetic 
field lines, a linkage inherited from the W-condensate linkages. 

This mechanism for magnetic field generation is somewhat different from 
Vachaspati's; it will be discussed in more detail in Section II. Although we 
do not even have the same scaling with the EW coupling g that Vachaspati 
does, we find similar numerical values for the magnetic field, of some 10 23 G 
at a magnetic scale length (transverse vortex size) Im ~ 2/g 2 T ~ 10~ 15 cm 
appropriate for the magnetic N = sector.^ 

The question of helicity generation is far less well-understood, because 
it is related to problems of formation and destruction of baryons (B) and 
leptons (L) in the early universe. We do not understand very well at all 
how B or L is generated in the early universe, but we do know something 
about how B+L can be destroyed by sphalerons and other EW effects at 
temperatures ranging from far above T c to just below it pT| , R|, [T3]. These 
EW effects all involve the B+L anomaly, so that changes in the number 
Nb,l are necessarily accompanied by changes in certain topological charges 

2 Numerical values are taken from Refs. Jl4|, 19, ^0|. It is natural, in view of the 
numerical values, to quote the magnetic mass and condensate values in terms of g 2 rather 
than the corresponding EW fine-structure constant aw = 3 2 /47r w 1/30. Other quantities 
are more naturally expressed in terms of the fine structure constants themselves, so a 
mixture of these appears in various expressions. 
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of the EW gauge theory. But further precision in this picture is hard to 
come by, since we do not presently know whether the present values of Nb 
and N L can be generated by EW effects, possibly including effects of bubbles; 
whether it is due to primordial (and conserved) B-L generation; what the net 
balance between creation and destruction by EW effects is, and so on. Nor 
do we understand in detail how helicity generated above T c can survive to 
below T c (but we will make some comments later in the spirit of the work of 



Martin and Davis |22| on the stability of magnetic fields generated at phase 
transitions). There are several mechanisms for translating the kind of EW 
helicity, or Chern-Simons terms, formed by B+L violation into magnetic 
helicity (of course, the usual Maxwell helicity does not occur in the B,L 
anomalies, since these currents are vectorlike with respect to the Maxwell 
field). We will discuss these mechanisms briefly in Section III, and postpone 
further elaboration to work now in preparation. For the present we will 
simply assume that the magnetic helicity surviving the EW phase transition 
is, give or take a couple of orders of magnitude: 

H M = a~\N B + N L ) ~10 66 erg- cm. (2) 

The numerical value comes from taking Nb ~ 10 80 ; a factor like a" 1 , where 
a is the fine-structure constant, is reasonable from the form of the B + L 
anomaly. 

Another way of thinking about Hm is to write its density as 

A-B = a~ 1 eT 3 (3) 

where the number density of any particle (or antiparticle) is ~ T 3 and e 
is a small number presumably related to, but not necessarily equal to, the 
B+L-entropy ratio (about 10~ 10 today[]23]]). In view of the smallness of e it 
may be wondered whether the helicity we assume can have any important 
effects. We argue that it does, on several grounds. First, up to about the 
time of decoupling helicity on large spatial scales is very nearly conserved, 
and at about this time the primordial helicity is likely to be large compared to 
natural scales for helicity in the universe then. In fact, although the helicity 
density is small compared to T 3 it is 30 or so orders of magnitude greater than 
one would expect from random fluctuations of the primordial magnetic fields 
in the absence of B + L- violating effects (just as the present B + L number 
is about 10 35 times the value expected from fluctuations alone). Second, the 
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presence of a primordial helicity can greatly affect the evolution of magnetic 
fields in the early universe, in particular by helping to generate an inverse 
cascade |T^, in which the original magnetic fields at the time of the 

EW phase transition, which have a very short magnetic correlation length 
l M ~ 10~ 15 cm are transformed into fields on much longer scales. (At the 
same time, there is a certain amount of destruction of the short-scale fields 
by magnetic viscosity; the universe is not a perfect conductor.) 

In Section IV we take up, in a simple approximation, the MHD issues 
concerning unstable dynamo growth driven by helicity in the early universe. 
This approximation is similar in spirit to two-scale models or mean-field 
models (see, e.g., Ref. m which the effect of fields fluctuating on small 

scales on long-scale fields is estimated by a process of spatially averaging 
quantities quadratic in the short-scale fields. Our approximation is such 
that magnetic helicity Hm is exactly conserved in the infinite conductivity 
approximation; finite-viscosity effects are trivially incorporated. We find a 
quasi-linear instability in which helicity transfers itself and the magnetic 
fields from short to long scales. There is an inverse- cascade instability for all 
lengths greater than a critical length scaling with e _1 , much larger than the 
magnetic length Im and much smaller than the Hubble size. It is perhaps 
surprising to find that although the maximum growth rate associated with 
the inverse cascade of helicity is 0(e 2 ), this rate can be comparable to the 
expansion rate of the universe. Part of this comes from the slowness of the 
expansion of the universe compared to the natural EW rates (EW rates are 
~ 10 14 times the expansion rate), and part from the occurrence of inverse 
powers of fine-structure constants in the formulas. 

Eventually the system approaches an equilibrium of the sort suggested in 
Refs. HTT| , |25| , in which the final scale length is ~ 1/e. This length is much 
bigger than Im but much smaller than the Hubble size. 

If the primordial helicity is large compared to fluctuation effects, it is 
interesting to ask whether it could be measured. We discuss this in Section 
V, and find it to be very doubtful, on two counts. First, measuring helicity 
of whatever size would require not only very complete Faraday rotation data, 
but also data on the angular gradients of the magnetic field; these would be 
very hard to get, but in principle are available from scattering from polarized 
dust grains. Second, the primordial helicity would dominate other effects up 
to about the time of structure formation, but after that time there are so 
many effects which process pre-existing magnetic fields that only very shaky 
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conclusions could be drawn. 



II. THE MAGNETIC EW CONDENSATE 



The general subject of the EW SU(2) condensate has been considered 
in detail elsewhere H [14]]. At temperatures above T c the Higgs field has 
vanishing VEV so the magnetic Ws are perturbatively massless in the N = 
Matsubara sector (in other Matsubara sectors the W-bosons either have an 
effective thermal mass ~ T or, in the N = electric sector a perturbative 
mass ~ gT) . However, it is well-known that non-perturbative effects generate 
a magnetic mass M w ~ g 2 T for the SU(2) gauge bosons, or Ws (but no 
such mass is generated for the hypercharge bosons). This mass is directly 
associated |TJ ; |26|] with a magnetic W-condensate 



(0) = \({G a l3 f) ~ g 6 T\ (4) 

and a negative free-energy density for the N = magnetic sector of value 
— (6)/ 3. Numerical estimates |TJ, |27|, |23|] for the free-energy density are in 
the neighborhood of —(0.01 — 0.02)g 6 T 4 . Because the neutral W potential is 
related to the usual Maxwell potential by A4 = sin OwWj* + cos 9wYi, (Yi is 
the hypercharge potential) there is a condensate of Maxwell fields too, with 
(including a factor of 4n for cgs units): 

(B 2 ) ~ 0.5e 2 /T 4 (5) 

and an RMS Maxwell field strength of about 10 23 G at T ~ T c . As already 
mentioned, the correlation length of this field above the EW phase transition 
is Im — 2/g 2 T ~ 10 -15 cm at T ~ T c . (For comparison, the EW Hubble 
scale is about 1 cm.) 

This field strength at the phase transition is numerically comparable to 
Vachaspati's estimate 0, but the scalings seem to be different. This is 
because he takes the W-mass to scale like gT instead of g 2 T. Both the present 
work and Vachaspati agree that B ~ Mw/g, but we use the usual magnetic- 
mass scaling for Mw. The numerical values agree fairly well because g ~ 0.65 
is not small compared to unity. 
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The spatial structure of this field is of interest. It could have finite-T 
sphalerons as a component, but as mentioned earlier [H| the mass of the 
magnetic sphaleron is so large that its Boltzmann factor is likely to be rather 
small, and we will ignore sphalerons. The dominant component is a gas 
of closed vortices of thickness M^ 1 whose entropy dominates the internal 
energy, so the vortices or strings like to be very long. They can link with 
each other and with themselves^, and also with Wilson loops; this latter 
is the mechanism for confinement both in three and four dimensions (in 
d = 4 the vortices are closed two-surfaces) p9[ , and the string tension is an 
outgrowth of fluctuations of the linking numbers of the vacuum condensate 
with the Wilson loop. The sum of all linking numbers of the vacuum (plus 
the sphaleron number, if any) is proportional to the Chern-Simons number of 
the thermal vacuum H [I4[ . This will be zero in the absence of such parity- 
violating effects as B+L violation; in its presence, the vacuum will have a 
net linking number. It is this linking number that we argue will be (in part) 
preserved as a contribution to Maxwell helicity after the EW phase transition 
is complete. 

We postpone a complete discussion of the transition dynamics to a later 
work. The general picture is that as the W-condensate dissolves (and the 
Higgs condensate forms), the phase- coherent Ws forming the condensate lose 
this coherence and become particles. We therefore need a picture of the 
condensate in particle language. Think first of a liquid-helium condensate: 
The atoms forming it are (in the center-of-momentum frame) strictly at zero 
momentum, even though the condensate coexists with other particles of finite 
momentum. The EW condensate is a little different because of the spatial 
structure of the vortices, which requires the Ws to have finite momentum, 
but smaller than the thermal momentum ~ T. In plasma language the vortex 
somewhat resembles a #-pinch, with pressure gradients balancing JxB forces. 
By looking at the vortex solutions (see the work of Cornwall in Ref. f29|j ) 
one reads off the scalings for SU (2) magnetic field B and condensate current 
density J and then infers the condensate pressure p: 

B~M 2 w /g; J ~ M^/g; p~g e T 4 . (6) 

The current is also ng, where n is the density of particles of one sign of charge, 
and the current velocity is of order unity. With Mw ~ g 2 T, one finds that 

3 Self-linking is also known as twisting and writhing; see for a discussion and 
references. 
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n ~ g 4 T 3 , so that a relatively small fraction of the particles in the thermal 
bath are tied up in the condensate. The particle energy in the condensate 
is nMw ~ g 6 T 4 , like the magnetic energy in Eq. (4). The Larmor radius 
of a particle scales like M w /gB ~ l/g 2 T = Im- The condensate particles 
and fields are in equilibrium above the phase transition, and depart from it 
relatively little during and just after the phase transition, in part because 
the W mass does not change very much during the transition. Of course, the 
photon mass does change drastically^ and new magnetic fields are generated, 
but their currents have for a while their original correlation length Im- 



III. B+L- VIOLATING GENERATION OF MAGNETIC 

HELICITY 



We are ignorant about the actual mechanism of B+L generation in the 
universe, although many candidates abound. The best we can do at present 
is to construct a scenario which is not obviously wrong; this scenario will not 
have precise numbers associated with the B+L violation effects. Imagine, 
then, that EW effects are well-described by the Standard Model, with no 
supersymmetry or other non-standard mechanisms being important. We 
will assume that the EW phase transition is second order (meaning a Higgs 
mass larger than about 80 GeV). In this case, there is no B+L production 
and no turbulence associated with the EW phase transition. The EW effects 
are limited to dissipation of some part of whatever B+L has been previously 
generated by some unknown mechanism. To avoid fine tuning, we assume 
that the original amount of B+L is comparable to what exists after the phase 
transition, or otherwise there is a cancellation between a large production rate 
and a large dissipation rate, leaving a much smaller net amount. By this 
assumption the amount of B+L which has been dissipated by EW effects is 
comparable to what exists today. 

We can roughly model these dissipation effects and their consequent gen- 
eration of topological charge by introducing a chemical potential associated 
with the conserved (but gauge- variant) current which is the sum of the B+L 

4 In fact, there is a sense (see Section 4) in which the photon mass becomes very slightly 
tachyonic, of O(e) . 
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current and a topological current: 



K» = J»+ ^^{Tr[g 2 A v d a A p - 9 -A u A a A p ] + \g a Y v d a Y^ (7) 
in standard notation. The grand partition function 

Z = Tre-^ H -" Q) , (8) 

where Q = J d 3 xK is the conserved charge, reinstates invariance under 
large gauge transformations by summing over all possible values of Nb+l- 
The chemical potential is of order n B+L T~ 3 , that is, 0(e). As far as the 
N = magnetic Matsubara sector of the SU{2) gauge fields is concerned, 
this chemical potential introduces a new gauge-coupling term ~ e which is 
CP-odd. 

There is, of course, no term in the B+L charge which directly involves the 
Maxwell helicity / d 3 xA ■ B, because of the vectorlike nature of electromeg- 
netism. But there are small terms in [iQ linearly coupled to A, which can 
be treated as perturbations on the underlying SU{2) condensate, in which 
electromagnetism participates. Since there is no hypercharge condensate^ 
we can just set the hypercharge potential to zero, so that one component 
(conventionally the third) of Wj is smOy/Ai. This is, at T > T c , a massive 
field with vortex solutions like the Nielsen- Olesen vortices; these solutions 
have been discussed in several places p9| |3|, 0, [31]. Let us write the effective 



CP-odd action for the Maxwell field as 



J d 3 xP ■ A. (9) 



Given that the underlying SU(2) condensate is one of strings, it turns out 
that P is effectively a magnetic field of the type 



Pi 



Y,qfd Zi M 2 A M (x-z) (10) 



where M (previously denoted M\y) is the magnetic mass, Aj^ is the massive 
Euclidean propagator for mass M, the integral runs over a closed string, the 
sum is over the collective coordinates of the strings in the condensate, and q 



5 At least, as driven by standard model fields; the Y boson remains massless at T > T c . 
Other sources of a Y condensate have been envisaged; see Joyce and Shaposhnikov |30[. 
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summarizes various constants and group matrices of no concern to us now. 
In a gauge where V • A = 0, the form of A to first order in e is: 



Mx) = c^^jkdj f dz k {A M - A }(x-z) + (V 2 + M 2 )" 1 ^ (11) 

(A is the massless propagator). Again the sum is over the strings of the con- 
densate, and c is a collection of constants of no interest now. This expression 



for A is reminiscent of a similar expression |y| for a vortex in d = 3 Yang- 
Mills theory with a Chern-Simons term added. Even if such a vortex string 
is straight one finds that it has a link number, coming from the twist of the 
field lines. Similarly, one readily calculates the magnetic helicity / d 3 xA ■ B 
from Eq. (11) to find that there is an 0(e) term from the product of the 
piece and the Pi piece, which, in its structure of string integrals, resembles a 
sum over regulated Gauss linking numbers. 

The calculation we have sketched out here relates the net SU(2) helicity, 
coming from B+L violation, to the currents driving the Maxwell magnetic 
field. It holds for T > T c . The next question to ask is how much of the 
magnetic helicity can survive the EW phase transition. If the problem is 
restricted to asking how much of the magnetic field, regardless of its helic- 



ity, survives, an answer has already been given by Martin and Davis [^2 
Their answer is that the field does survive largely intact. Although we will 
not discuss it in detail, we believe that most of the helicity also survives. 
One approach is similar to that of Martin and Davis: Just below the phase 
transition, when the Higgs field (f) has a VEV (denoted v) one can use M 



't Hooft's [32 expression for the Maxwell field strength involving projecting 



out the Abelian part with the Higgs field. This can be written: 

Ai 

F l3 = diAj - — sin 9w(di<l>)%<j> - (i <-> j). (12) 

where the electromagnetic vector potential is defined in terms of the unit 
Higgs vector as: 

Ai = sia6 w n a W % a ; n a = ^a a (j)/v (13) 

(There is also an electric field which we will not consider here.) The phase of 
the Higgs field near the transition temperature inherits the topological infor- 
mation on the linkage of the string condensate. Martin and Davis consider 
thermal fluctuations at a second-order phase transition and conclude that 
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these do not wipe out the magnetic field. The same conclusion ought to hold 
for the helicity as well. 

Another purely classical argument has been given by Taylor \1E\ for Toka- 
mak plasmas, which have macroscopically linked magnetic fields. He begins 
with a non-equilibrium plasma, of the sort that might be produced at a phase 
transition, and points out that for closed field lines in the ideal MHD case, 
the magnetic helicity ought to be preserved separately for every closed field 
line (or surface). He then argues that as the plasma relaxes to equilibrium 
through effects due to large but finite conductivity and field lines reconnect, 
the helicities on the various field lines will be homogenized, but the sum (i.e., 
volume integral) of all helicities should be preserved, because reconnection — 
which homogenizes helicity — does not make large changes in the fields and 
potentials themselves. 

Given that some helicity of order e survives, how does it affect the evolu- 
tion of magnetic fields after the EW transition? We discuss this next. 



IV. MAGNETOHYDRO DYNAMICS AND MAGNETIC 

HELICITY 



Although it is known [y, [10|, [TJ, |1J that helicity is an important driver 
of unstable MHD dynamo action, the problem we face is somewhat different 
from the standard MHD dynamo problem because general relativity must 
be taken into account |24]] and because we consider an initial-value problem 
rather than the often-considered problems of steadily-driven helicity. 

Let us consider some typical length scales, beginning with the Larmor 
radii. In the SU(2) condensate above the phase transition, the coherent 
Larmor radius Rl, or the radius a charged particle in the condensate would 
have if circulating around a single vortex, is about 

T>T C : R L ps g 2 T/gB ~ \jg 2 T\ (14) 

a non-condensate particle would have a Larmor radius bigger by only a few, 
a factor of 1/g 2 , modified possibly for the fact that it is circulating around 
a couple of vortices so the RMS B-field is a little smaller than the coherent 
field. In the plasma after the phase transition some factors of g change to 
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factors of e (compare Eqs. (4) and (5)), and the former condensate particles 
have a Larmor radius 

T<T C : R L &l/e 2 T, (15) 

a few times larger than before the phase transition. These Larmor radii are 
bigger than Coulomb collision lengths, which scale like l/a 2 T (for quarks the 
QCD length is about l/a|T; these could be unmagnetized, but they play no 
role in the condensate formation anyhow). In what follows, we will assume 
that the Larmor radius of Eq. (15) is appropriate for the inital scale lengths 
of Maxwell fields just after the phase transition. 

The helicity-driven dynamo picture is that these small-scale fields drive 



an instability of fields on larger scales. As is usual in such cases |]IT| , [12], [13 
we will average quantities quadratic in fields varying on small scales to get 
driving terms on the large scales. 

After the modifications for general relativity our approach is similar to 



that of Ref. [TT|] for so-called a dynamos, but with some significant differ- 
ences. Brandenburg, Enqvist, and Olesen [|4j] have shown that the ideal 
MHD equations (no kinetic or magnetic viscosity) have simple scaling prop- 
erties in a flat Robertson- Walker metric 

ds 2 = dt 2 - R 2 {t)dx 2 (16) 

which reduces them to the flat-space MHD equations provided that cosmic 
time t is replaced by conformal time t: 

i= J dt/R(t) = 2{t EW t) 1/2 (17) 

where we normalize so that at EW time tEw ~ 10~ n sec the scale factor 
R(t) is unity. The required ideal-MHD scalings are: 

B = B/R 2 } J = J/R 3 , p = p/R\ v = v. (18) 

Here B,J,p,v are respectively the magnetic field, the current, the MHD 
energy density, and the bulk velocity. Other scalings can be inferred from 
these and the statement that only time is changed, not spatial gradients. The 
meaning of these scalings is that the variables with a tilde do not change just 
because the scale factor R is changing. The scaling of magnetic field reflects 
the usual dilution with flux conserved. 
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We need two additional scalings, since we will consider non-ideal MHD. 
These scalings are for the conductivity a and the collision (or correlation) 
time r. In the case of Coulomb collisions it is well-known that 

a « T/a = a/R; r c « (a 2 T)- 1 = f c R. (19) 

The magnetic viscosity vu is 1/a and scales like R; so does the kinetic 
viscosity (and also a kind of Alfven viscosity defined later). 

If one writes out the general-relativistic MHD equations in terms of the 
original variables one finds that they are the same as the flat-space MHD 
equations with the tilde variables. Before we write these equations down we 
make two points: First, the bulk velocity v is expected to be non-relativistic, 
even if the phase transition is first-order. Of course, the thermal velocities 
of the (non-condensate) plasma are essentially the speed of light, so the 
energy density p is quite relativistic. Second, we will assume a second-order 
transition, which means that there should be no turbulence in the matter 
variables and that, because of inflation, the scale lengths for these variables 
are quite large, at least initially.^ As the instability we identify below grows, 
the matter variables at later times will also become turbulent, which can 
have important effects |13| not considered here. 

To avoid notational complexity we drop the tildes on the scaled variables 
except for that denoting conformal time. We also drop all gradients referring 
to matter variables, including the pressure and kinematic viscosity terms as 
well as the non-linear velocity term.[] For the moment we will keep magnetic 
viscosity. Then the matter equation is (using J = V x B) 

(4/3)p^ = -B x (V x B). (20) 

Using the electric-field relation 

E = -vx B + J/a (21) 

6 But the scale lengths for the current J are the same as for the magnetic field and 
helicity, that is, of 0(1/ e 2 T). There is no reason for the scale length for other matter 
variables to be this short. 

7 It is only the pressure gradient that is dropped; the pressure contributes to the overall 
energy density. 
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we can strip off a curl operation from the B equation to get: 

8 A 1 - - 

— = v x B - -V x B + Vy. (22) 
dt a K ' 

The Vx term has no physical effect and can be dropped. 

Magnetic viscosity is due to the 1/a term, and is dominated by Coulomb 
scattering (particle viscosity is dominated by quark scattering). A numerical 
estimate at the EW phase transition gives 

v M = 1/a « 1(T 8 - 10~ 9 cm 2 /sec. (23) 

The corresponding magnetic Reynolds number Iv/vm ranges from about 10 5 
at the correlation length scale I ~ 1/ay/T to 10 19 at the Hubble scale I = 
H" 1 , assuming that the velocity is the speed of light (actually it is somewhat 
smaller). We will identify below a sort of Alfven magnetic viscosity which 
is somewhat larger than the collisional viscosity, and therefore drop the 1/a 
term in Eq. (22). 

The next step in principle is to solve for the velocity from Eq. (20). Of 
course this cannot be done analytically, and we make a common approxi- 
mation (see, e.g., Ref. |L3], |33|]) and invert the i integration in Eq. (20) 
by multiplying by a correlation time r. It is reasonable to choose r to be 
the Coulomb time r c w l/a 2 T (most particles do not change their veloc- 
ity appreciably in one Larmor period). Assuming that the initial velocity 
vanishes: 

v = ~B x (V x B). (24) 

— * 

Then the equation for A is: 

°-~ = 6 ^{B[B ■ (V x B)) - (V x B)B*}. (25) 

This highly-nonlinear equation is intractable as it stands, so we quasi- 
linearize it by replacing two of the three factors involving the magnetic field 
on right-hand side by averages (in the spirit of the mean-field dynamo [|TT 



in 



which averages are made over short spatial scales to find their effect on long 
scales). We assume, by isotropy, that all vectorial quantities, in particular 
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B, have zero average, so the only possibility is to replace Eq. (25) by: 



f)A -, -> -, -> -> -> 

-g = j^{B(B • (V x B)) - (V x B)(B 2 )}. (26) 

This is of the canonical form of an a-dynamo with turbulent viscosity [ pi]] : 

dA -> - 

- r = a M B-(3 M (VxB). (27) 
at 

So far we have not been specific about the nature of these averages. To 
do that, we note that once having dropped the magnetic viscosity, it should 
be that the magnetic helicity is conserved. Compute the rate of change of 
helicity in a large volume V: 

dA 

H M = 2 / d 3 x—^ ■ B 



v dt 



= 2 f d 3 x{a M B 2 - p M [B ■ (V x B)]}. (28) 
Clearly if one defines averages so that 

a ^ = ^r\ d " x ^ ■ (V x B), p M = [ d 3 xB 2 (29) 
41/ p Jv AV p Jv 

then helicity is conserved. 

One may now go back and add the collisional magnetic viscosity vm to 
(3m m Eq. (27), an addition which will lead to destruction of helicity on short 
scales. We will not consider the collisional magnetic viscosity further. 

Although ckm is not the magnetic helicity Hm = I A ■ B defined earlier, 
it is qualitatively equivalent to it, for the initial magnetic field configuration 
from which we start, which is dominated by a single scale length Im or Rl. 
We will estimate the coefficients qm, (3m just at the EW transition point. 
At later times these values will change, not only because of i?-scaling but 
also because of MHD effects. Let us estimate <xm m terms of H M - To find 
(B-(VxB)) we divide H m by the square of the Larmor length Rl 1/ e 2 T. 
This yields: 

r A ■ B ee 4 
a M ~ - — =— ~ 0.01— w 3 x 10 12 e cm/sec. (30) 
pR L a 6 
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(3 M ~ ^(B 2 ) ~ U ' « 10- 6 cm 2 /sec, (31) 



(The numerical factors in Eqs. (30) and (31) below summarize various factors 
involved in the densities, etc; we are using units such that e ~ 0.3, g ~ 0.65.) 
For p M : 

r ^ 0« 
p cri 

considerably larger than the collisional magnetic viscosity of Eq. (23). The 
ratio Pm/c^m defines a critical length, which we call In, separating the larger 
scales to which helicity is transferred from the smaller scale from which it 
comes. This is approximately 3 x 10~ 18 e -1 cm. So In is much larger than 
any correlation length but much smaller than the Hubble size. 

Note that (3m can be written as t ( V\ in terms of a (nominal) Alfven 
velocity Va ~ eg 2 . This suggests that the /^-effect can be interpreted in 
terms of Alfven waves carrying off helicity from small scales to large. One 
can, in fact, easily check that Alfven waves do carry helicity, but it is not 
clear that the Alfven wave concept is suited to a plasma of the type found 
in the early universe. 

Next we study the quasi-linear instability resulting from equation Eq. 
(27), replacing qm, Pm by constants. In fact, these quantities change, so the 
approximation of constant values can only hold for a short time. There have 
been numerous studies of a-dynamos, described by equations like Eq. (27). 
We give a treatment slightly different from any of which we know. Begin 
with the standard [11| decomposition of the magnetic field into poloidal (P) 



and toroidal (T) components: 

B = LT + V x (LP) (32) 
with the magnetic potential in a natural gauge: 

A = -rT + LP. (33) 

In these equations, 

L = fxV (34) 

is the anti-Hermitean generator of angular momentum. The helicity is given 
by: 

E M = I d?xA ■ B = -2 f d 3 xTL 2 P (35) 



and expresses the linkage between toroidal and poloidal field lines, both of 
which must be present for there to be helicity. (This fact is related to the 
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possibility of measuring cosmic helicity, which we discuss in Sec. V.) The 
equations for T, P are easily derived: 



f = -a M V 2 P + (3 M V 2 T (36) 

P = a M T + (3 M V 2 P (37) 

These are straightforwardly solved by Fourier transforms (indicated by a 
hat): 

T(x, t) = ^-y 3 J d 3 k[f + e x ^ + Tie*-*] (38) 

p ^>w/t^-^ (39) 

where 

f ± = f±P\ f=0 , \± = ±ka M -k 2 (3 M - (40) 

Whatever the sign of olm there is growth (unless the corresponding T iden- 
tically vanishes); to be definite assume a« > 0. The maximum growth rate 
occurs for k ~ Ijj 1 , and is given by 

2 

A = Max A + = ^^ 3 x 10 30 e 2 sec" 1 . (41) 

4/3 M 

Since the expansion rate of the universe H is about 10 11 sec -1 , there will be 
appreciable dynamo growth for e > 2 x 10~ 10 , a not unreasonable number.^ 
As Ref. |Tl[] has pointed out, there are also equilibrium solutions to the 
equations (36, 37). Naturally, their scale length is about In- The equilibrium 
solution for P or T is tachyonic: 

V 2 P+(l H y 2 P = 0. (42) 

Taylor [^] has given another understanding of this tachyonic equilibrium 
equation, by introducing a real Lagrange multiplier corresponding to the 
conserved helicity Hm, which is a Chern-Simons term with real coefficient. 
It is well-known that such a system is tachyonic. However, no real problems 
are encountered; the solutions to Eq. (42), involving Bessel functions (for 

8 Recall that just after EW times the fractional baryon number is about 3 x 10 -8 , rather 
larger than it is now, because more photons have been produced at annihilation |23| . 
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a given angular momentum), can be terminated smoothly matching on to 



multipolar magnetic fields in the region where the helicity vanishes [J I 



It is not possible for us to go to large times (say, larger than the expansion 
time of the universe) by analytical means, because the coefficients «m, Pm 
are changing as the magnetic scale lengths change, and as magnetic energy 
is exchanged with fluid energy. We hope to carry out elsewhere a detailed 
numerical simulation of the MHD processes so crudely described here. This 
simulation, going to later times, would reveal the influence of a number of 
effects we have omitted, including the growth of kinetic helicity in the MHD 
velocity field, as measured by J pv ■ (V x v). This should grow from its 
initial value of zero to approximate equipartition with the magnetic helicity 
as measured by / B ■ (V x B), at which point one encounters the so-called 
Alfven effect of Ref. fl3f : These two helicities tend to cancel out on short 
spatial scales, leaving still a residual inverse-cascade instability. 



V. CAN WE DETERMINE THE EXISTENCE OF 
PRIMORDIAL HELICITY EXPERIMENTALLY? 



We remark on our assumed scale for magnetic helicity (see (2)) in com- 
parison with the value expected from fluctuations at EW time, and on the 
possibility of measuring helicity today. 

Supposes that at EW time the net helicity of the universe were given 
only by fluctuations of random-sign helicity in cubes the size of the magnetic 
correlation length Im- The maximum helicity in any cube scales like (B 2 )l 4 Ml 
so the RMS fluctuation value is found by multiplying by the square root of 
N « V/l\ 4 , where V is the volume of the universe at EW time. Taking (B 2 ) 
from (5) one finds, with V ~ 10 42 cm 3 , that the helicity due to fluctuations is 
about 10 30 erg-cm, far smaller than we estimated in equation (2). This is, of 
course, about the same as for the actual number of baryons today compared 
to the fluctuation value of about 10 45 . 

If we could measure the EW-time helicity it would be easily possible 
to prove or disprove our hypothesis that there is a primordial helicity pro- 
portional to the baryon number. Unfortunately we cannot, and the enor- 
mous processing of magnetic fields that must have taken place after matter 
dominance and structure formation would in any case greatly obscure any 
interpretation of a measurement of helicity today (that is, since structure 
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formation). But it may still make some sense to ask whether we could in 
principle measure the present magnetic helicity of the universe. 

We have seen from Eq. (35) that to measure helicity requires simulta- 
neous determination of both T and P in the decomposition Eq. (32) of the 
magnetic field. Faraday rotation can, in principle at least, recover the com- 
ponent r ■ B, equivalent to P (given enough sight lines to polarized sources 
and independent measurements of the electron density). However, a sepa- 

— * — * 

rate measurement is needed, which must amount to the extraction of L ■ B, 
equivalent to T. Again in principle, this could be measured from polarization 
of starlight by scattering from dust grains polarized in the cosmic magnetic 
field, given enough separate dust clouds, and a hypothesis of general isotropy 
of the magnetic field. 

Finally, one might think of measuring a net circular polarization in the 
cosmic microwave background, but if this is only 0(e) this too is probably 
impossible. 

Even if we could do this, would the primordial helicity of Eq. (2) be 
large compared to fluctuation-driven helicity today, as it was at EW times? 
The answer is no, because even for a single galaxy with B ~ 1CT 6 G, an in- 
plane correlation scale of 15 kpc and a vertical correlation scale of 0.5 kpc 
the maximum helicity is ~ 10 77 erg-cm, far greater than a primordial value. 
This very large value on the primordial scale is analogous to the very large 
value of the fields and coherence lengths compared to those of primordial 
fields (scaled to the time of structure formation) with or without helicity, 
and simply reflects the well-known fact that today's fields grew exponen- 
tially during the time since structure formation, by dynamo action. (It is, 
of course, also possible that primordial fields and helicity have nothing to do 
with today's cosmic magnetic fields.), Still, measuring the helicity of today is 
interesting in principle, since it might very well be larger than one would ex- 
pect from fluctuations alone, even though there has been so much processing 
of magnetic fields since the beginning of structure formation. 
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